Weak localisation magnetoresistance and valley symmetry in graphene 
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Due to the chiral nature of electrons in a monolayer of graphite (graphene) one can expect weak 
antilocalisation and a positive weak-field magnetoresistance in it. However, trigonal warping (which 
breaks p -^ — p symmetry of the Fermi line in each valley) suppresses antilocalisation, while inter- 
valley scattering due to atomically sharp scatterers in a realistic graphene sheet or by edges in a 
narrow wire tends to restore conventional negative magnetoresistance. We show this by evaluating 
the dependence of the magnetoresistance of graphene on relaxation rates associated with various 
possible ways of breaking a 'hidden' valley symmetry of the system. 
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PACS numbers: 73.63.Bd, 71.70.Di, 73.43.Cd, 81.05. Uv 

The chiral nature 0, S S Hi of quasiparticles in 
graphene (monolayer of graphite), which originates from 
its honeycomb lattice structure and is revealed in quan- 
tum Hall effect measurements ^ ^, is attracting a lot 
of interest. In recently developed graphene-based tran- 
sistors la, |a| the electronic Fermi line consists of two 
tiny circles Q surrounding corners K.± of the hexagonal 
Brillouin zone Q, and quasiparticles are described by 4- 
component Bloch functions <& =[4>-k..a^ 4>-k..bj 'Pk^ b^ 
<Pi^_ ^], which characterise electronic amplitudes on two 
crystalline sublattices {A and B), and the Hamiltonian 

H ^vllz®crp- fi[(7x{pl --py) - 2ayPxPy]. (1) 

Here, we use direct products of Pauli matrices 
fa;,j,,z, ctq = 1 acting in the sublattice space {A, B) and 
nx,j/,z,no = i acting in the valley space (K±) to high- 
light the form of H in the non-equivalent valleys [g • Near 
the center of each valley electron dispersion is determined 
by the Dirac-type part v crp of H. It is isotropic and lin- 
ear. For the valley K+ the electronic excitations with mo- 
mentum p have energy vp and are chiral with crp/p — 1, 
while for holes the energy is —vp and crp/p = — 1. In the 
valley K_, the chirality is inverted: it is crp/p = — 1 for 
electrons and crp/p = 1 for holes. The quadratic term in 
Eq. 1^ violates the isotropy of the Dirac spectrum and 
causes a weak trigonal warping 8] . 

Due to the chirality of electrons in a graphene-based 
transistor, charges trapped in the substrate or on its sur- 
face cannot scatter carriers in exactly the backwards di- 
rection I2, Qj provided that they are remote from the 
graphene sheet by more than the lattice constant. In 
the theory of quantum transport Q the suppression of 
backscattering is associated with weak anti-localisation 
(WAL) f^. For purely potential scattering, possible 
WAL in graphene has recently been related to the Berry 
phase TT specific to the Dirac fermions, though it has also 
been noticed that conventional weak localisation (WL) 
may be restored by intervalley scattering [l^, uM ■ 

In this Letter we show that the WL magnetoresistance 



in graphene directly reflects the degree of valley sym- 
metry breaking by the warping term in the free-electron 
Hamiltonian ^ and by atomically sharp disorder. To 
describe the valley symmetry, we introduce two sets of 
4x4 Hermitian matrices: 'isospin' E = {T^x,^y,^z) and 
'pseudospin' A = (A^;, Ay, A^). These are defined as 
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= H^ <Siao, 


(3) 



and form two mutually independent algebras, [E, A] = 0, 

[E,,,E,J=2i£^i^^^S„ [Ai„Ai,]^2ie'''^'Ai, 

which determine two commuting subgroups of the group 
U4 of unitary transformations [l3| of a 4-componcnt $: 
an isospin (sublattice) group SUl" = {e*°"'^} and a pseu- 
dospin (valley) group SU^ = {e**"'^}. 

The operators E and A help us to represent the electron 
Hamiltonian in weakly disordered graphene as 

H = v^p + h„ + iu{r)+ Y^ EsA,u,,i(r), (4) 

s,l—x,y,z 



where h^ 



-/iE,(Ep)A,E,(Ep)E, 



The Dirac part of H in Eq. I^J , v Ep and potential dis- 
order Iu(r) [I is a 4x4 unit matrix and {u{r)u{r')) — 
v?5 (r — r')] do not contain pseudospin operators A/, i.e., 
they remain invariant under the group SU2 transfor- 
mations. Since E and A change sign under the time- 
inversion |l4| . the products E^A/ are t — > — i invariant 
and, together with I can be used as a basis to repre- 
sent non-magnetic static disorder. Below, we assume 
that remote charges dominate the elastic scattering rate, 
T^^ f^ Tq = Ti^v? /h, where 7 = p-p / [2^11? v) is the 
density of states of quasiparticles per spin in one valley. 
All other types of disorder which originate from atom- 
ically sharp defects [13 and break the SU^ pseudospin 
symmetry are included in a time- inversion-symmetric |l4| 



random matrix T,sAiUsj{r). Here, ^^.^(r) describes dif- 
ferent on-site energies on the A and B sublattices. Terms 
with Ua:,z(r) and Mj,_z(r) take into account fluctuations of 
A ^ B hopping, whereas Us,a;(r) and Us^y(r) generate 
inter-valley scattering. In addition, warping term h^, not 
only breaks p — > — p symmetry of the Fermi lines within 
each valley but also partially lifts SU^-symmetry. 

Hidden SU2 symmetry of the dominant part of H in 
Eq. I0J enables us to classify the two-particle correla- 
tion functions, 'Cooperons' which determine the inter- 
ference correction to the conductivity, 6g by pseudospin. 
Below, we show that Sg is determined by the interplay 
of one pseudospin singlet (C^) and three triplet (C^'^^^) 
Cooperons, Sg oc — C° + C^ + C^ + C^ , some of which 
are suppressed due to a lower symmetry of the Hamilto- 
nian in real graphene structures. That is, the 'warping' 
term /iw and the disorder E^A^Us^^ suppress intravalley 
Cooperons C^'^ and wash out the Berry phase effect and 
WAL, whereas intervalley disorder SgAj-jj^jMs ,j,(j^)(r) sup- 
presses C^ and restores weak localisation 9] of electrons, 
provided that their phase coherence is long. This results 
in a WL-type negative weak field magnetoresistance in 
graphene, which is absent when the intervalley scatter- 
ing time is long, as we discuss at the end of this Letter. 

To describe quantum transport of 2D electrons in 
graphene we (a) evaluate the disorder-averaged one- 
particle Green functions, vertex corrections, Drude con- 
ductivity and transport time; (b) classify Cooperon 
modes and derive equations for those which are gapless in 
the limit of purely potential disorder; (c) analyse 'Hikami 
boxes' [^ [I3 for the weak localisation diagrams paying 
attention to a peculiar form of the current operator for 
Dirac electrons and evalute the interference correction to 
conductivity leading to the WL magnetoresistance. In 
these calculations, we treat trigonal warping h„ in the 
free-electron Hamiltonian Eqs. (|1I4I) perturbatively, as- 
sume that potential disorder lu{r) dominates in the elas- 
tic scattering rate, t~^ '^ Tq = tt'ju^ /h, and take into 
account all other types of disorder when we determine 
the relaxation spectra of low-gap Cooperons. 

(a) . Standard methods of the diagrammatic technique 
for disordered systems J|, TQj at ppVT ^ h yield the 
disorder averaged single particle Green's function, 



G-M(P,.) = ^|^^±^ 



4yA-^V' '"^^^ 



= e±^ihTf^\ 



The current operator, v = vT. for the Dirac-type par- 
ticles described in Eq. |^ is a momentum- independent. 
As a result, the current vertex Vj { j — x,y), which enters 
the Drude conductivity. Fig. 1(a), 
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is/li^^-l*''^''*'''*''*''--"} 



= 4e^7D, with D = i/ro = ^v^tu, (5) 

is renormalised by vertex corrections in Fig. ^b) : v = 
2v = 2uS. Here 'Tr' stands for the trace over the AB 
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FIG. 1: (a) Diagram for the Drude conductivity with (b) 
the vertex correction, (c) Bethe-Salpeter equation for the 
Cooperon propagator with valley indices CmCV 8.nd AB lat- 
tice indices ajSa'P' . (d) Bare 'Hikami box' relating the con- 
ductivity correction to the Cooperon propagator with (e) and 
(f) dressed 'Hikami boxes'. Solid lines represent disorder av- 
eraged G^' , dashed lines represent disorder. 



and valley indices. The transport time in graphene is 
twice the scatering time, Ttr = 2tq, due to the scatter- 
ing anisotropy (lack of bacskattering off a potential scat- 
terer). This follows from the Einstein relation Eq. Q 
(where spin degeneracy has been taken into account). 

(b) . The WL correction to the conductivity is associ- 
ated with the disorder-averaged two-particle correlation 
function C'^g'^/^/ known as the Cooperon. It obeys the 
Bethe-Salpeter equation represented diagrammatically in 
Fig. njc). The shaded blocks in Fig.^c) are infinite se- 
ries of ladder diagrams, while the dashed lines represent 
the correlator of the disorder in Eq. Q). Here, the valley 
indices (K±) of the Dirac-type electron are included as 
superscripts with incoming ^/i and outgoing ^ /i', and the 
sublattice (AB) indices as subscripts a(3 and a'/?'. 

It is convenient to classify Cooperons in graphene as 
iso- and pseudospin singlets and triplets, 



1/3 



><C%%{^s.^yAiAy)';'i'- (6) 

Such a classification of modes is permitted by the com- 
mutation of the iso- and pseudospin operators S and 
A in Eqs. (|2EEIl, [^s,M = 0. To select the isospin 
singlet (s = 0) and triplet (s — x, y, z) Cooperon com- 
ponents (scalar and vector representation of the group 

{e""'^}), we project the incoming and outgo- 



SU? - r^"n-s 



ing Cooperon indices onto matrices SySs^and Hg^Yiy, re- 
spectively. The pseudospin singlet (/ = 0) and triplet 
(Z = X, y, z) Cooperons (scalar and vector representation 



of the 'valley' group SU^ = {e*^"'^}) are determined by 
the projection of C^^'^/i, onto matrices AyAj^ (Ai^Ay) 
and are accounted for by superscript indices in Cj^'j^ . 
For disorder Iu(r), the equation in Fig. QJc) is 



e(Kzp,p) [g. It has been noticed [T^ that such a de- 
formation of a Fermi line of 2D electrons suppresses 
Cooperons. As h^ has a similar effect, it suppresses the 
pseudospin-triplet intravalley components Cg and Cq , at 
the rate 
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obtaining the diffusion equations for the Cooperons using 
the gradient expansion of the Bethe-Salpeter equation we 
take into account its matrix structure. The matrix equa- 
tion for each set of four Cooperons C', where I = Q,x,y, z 
has the form 
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After the isospin-triplet modes were eliminated, the diffu- 
sion operator for each of the four gapless/low-gap modes 
Cq becomes Dq"^ — iw -I- Fg, where D = ^v'^Tti = v'^tq. 

Symmetry-breaking perturbations lead to relaxation 
gaps Fq in the otherwise gapless pseudospin-triplet com- 
ponents, Cq,Cq,Cq of the isospin-singlet Cooperon 
Cq, though they do not generate a relaxation of the 
pseudospin-singlet Cq protected by the time-reversal 
symmetry of the Hamiltonian Q . We include all scatter- 
ing mechanisms described in Eq. Q) in the correspond- 
ing disorder correlator (dashed line) on the r.h.s. of the 
Bethe-Salpeter equation and in the scattering rate in the 



;i = 2to {e^fi/hv^y 



(8) 



X Tr{l],I],AiA, [g^,,^+J A,A,,S,S,,G^q_p^,}. and Q 



However, since warping has an opposite effect on valleys 
K+ and K_, it does not cause gaps in the intervalley 
Cooperons Cq (the only true gapless Cooperon mode) 



It leads to a series of coupled equations for the Cooperon 
matrix C' with components C"^, . It turn out that for 
potential disorder lu{r) isospin-singlet modes Cgg arc 
gapless in all (singlet and triplet) pseudospin channels, 
whereas triplet modes C"^ and Cyy have relaxation gaps 

ri = Ty = ^Tq^ and Ci^^ have gaps F^ = Tq\ When ^ead 



Altogether, the relaxation of modes Cg can be de- 
scribed by the following combinations of rates: 

F° — n F^ — 9t^^ px _ py _ -1 I -1 I -1 ^ -1 



In the presence of an external magnetic field, B = rotA 
and inelastic dccoherence, r~^, equations for Cg = Cgg 



[D{^V + ^Af + F[, + r-i - ^c.]&Q (r, r') ^ S (r - r') . 

(c). Due to the momentum-independent form of the 
current operator v =2vT,, the WL correction to conduc- 
tivity Sg includes two additional diagrams, Fig.^^e) and 
(f) besides the standard diagram shown in Fig. QJd). 
Each of the diagrams in Fig. Qfe) and (f) [not included 
in the analysis in R.ef.jlll]] produces a contribution equal 
to (— -j) of that in Fig. ^d). This partial cancellation, 
together with a factor of four from the vertex corrections 
and a factor of two from spin degeneracy leads to 
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(9) 



Using Eq. ©, we find the B = temperature de- 
pendent correction, 5p to the graphene sheet resistance, 



5p{Q) e^ 

— ~ ^ ^Sg = — 

/9^ TTAl 



ln(l + 2^)-21n^^ 



, (10) 



disorder-averaged G^/^, as Tg ^ ^ t ^ =Tq^+Y.sI t si' ■ and evaluate magnetoresistance, p{B) - p{Q) = Ap{B), 



For simplicity, we assume that different types of disorder 
are uncorrelated, {us,i{r)us',i'{r')) = uliSss'Su'S{r - r') 
and, on average, isotropic in the x — y plane: u^j — u^i = 
= uj. = uj.^. We parametrize them by scatter- 
g rates t^;^ = TT^ulJh, where t^^ — rjy = tJI^ and 
/ ^; = T^i = t]]; due to the x — y plane isotropy of dis- 
order, which are combined into the intervalley scattering 
rate rr and the intra- valley rate t~^, as 



ApiT,B) 



eV 
nh 



F{^)-Fi 



B 



"■±1^ 
ini 



'B^ 



-2F{ 



Bin + 2B[ 



B 



B^ + _B* 



(11) 



F{z)^lnz + iPil + -), B^.i,, 
2 z 



he 

llTe 



4r; 



^r-l r-' 



4t^ 



2r-i. 



(7) 



The trigonal warping term, h^ in the Hamiltonian (^ 
plays a crucial role for the interference effects since it 
breaks the p — > — p symmetry of the Fermi lines within 
each valley: e(K±, — p) ^ e(K±,p), while e(K±, — p) = 



Here, ^ is the digamma function, and the decoherence 
tZ,^{T) determines the MR curvature aX B < B^. 

Equations (|lll) and H1U|I represent the main result of 
this paper. They show that in graphene samples with 
the intervalley time shorter than the decoherence time, 
Tip > Ti, the quantum correction to the conductivity has 
the WL sign. Such behavior is expected in graphene 
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FIG. 2: MR expected in a phase-coherent graphene t^ S> ti: 
with TzjTw ^ Ti (dashed) and r, ^ ri (solid hne). In the 
case of Tip < Ti, Sp — 0, so that Ap{B) = 0. 



tightly coupled to the substrate (which generates atom- 
ically sharp scatterers) . Figure |21 illustrates the corre- 
sponding MR in two regimes: B* ^ Bi [tz^t^ :s> Ti) and 
-B* ^ Bi (r, ^ Ti). In both cases, the low-field MR, at 
B < Si is negative (for B, ^ Bi, the MR changes sign at 
B ^ Bi). A dashed line shows what one would get upon 
neglecting the effect of warping, the solid curve shows 
the MR behavior in graphene with a high carrier density, 
where the effect of warping is strong and leads to a fast 
relaxation of intravalley Cooperons, at the rate described 
in Eq. (0). Then, in Eqs. H10lll|l r, w r^ < Ti < r^ and 
B* ^ Bi, which determines MR of a distinctly WL type. 
Note that in the latter case MR is saturated at i? ~ i?i, in 
contrast to the WL MR in conventional electron systems, 
where the logarithmic field dependence extends into the 
field range of hc/ADerti-. In a sheet loosely attached to a 
substrate (or suspended), the intervalley scattering time 
may be longer than the decoherence time, Ti > r^ > Tw 
(Bi < B^p < B^). In this case, Cq in Eq. © is effec- 
tively gapless and cancels Cq, whereas trigonal warping 
suppresses the modes Cq and Cg , so that 5g = Q and MR 
displays neither WL nor WAL behavior: Ap{B) = 0. 



Equation (|11() explains why in the recent experiments 
on the quantum transport in graphene |l8l | the observed 
low-field MR displayed a suppressed WL behavior rather 
than WAL. For all electron densities in the samples stud- 
ied in 18i] the estimated warping-induced relaxation time 
is rather short, T^/rtr ~ 5-^30, r^ < t^ which excluded 
any WAL. Moreover, the observation 1»| of a suppressed 
WL MR in devices with a tighter coupling to the sub- 
strate agrees with the behaviour expected in the case of 
sufficient intervalley scattering, Ti < t^, whereas the ab- 
sence of any WL MR, Ap{B) = for a loosely coupled 
graphene sheet is what we predict for samples with a long 
intervalley scattering time, Ti > t^. 

In a narrow wire with the transverse diffusion time 
L\/D <^ Ti,T*,T;p, edges scatter between valleys [l7| . 
Thus, we estimate Fq ~ tt'^D/L']_ for the pseudospin 
triplet in a wire, whereas the singlet Cq remains gapless. 
This yields negative magnetoresistivity for B < 2ttB±, 
B±_ = hc/eL\: 



APwirc {B) 2e2L, 



P' 



h 



1 + iB^B^Bj 



(12) 



Equations H1(JI12|I completely describe the WL effect in 
graphene and explain how the WL magnetoresistance re- 
flects the degree of valley symmetry breaking. They show 
that, despite the chiral nature of electrons in graphene 
suggestive of antilocalisation, their long-range propaga- 
tion in a real disordered material or a narrow wire does 
not manifest the chirality. 
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